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It is well-known that in bulk, the solution of the Bogoliubov-de Gennes equations is the same 
whether or not the Hartree-Fock term is included. In this case the Hartree-Fock potential is position 
independent and, so, gives the same contribution to both the single-electron energies and the Fermi 
level (the chemical potential). Thus, the single-electron energy measured from the Fermi level 
(it controls the solution) stays the same. It is not the case for nanostructured superconductors, 
where quantum confinement breaks the translational symmetry and results in a position dependent 
Hartree-Fock potential. Now the contribution of the Hartree-Fock mean field to the single-electron 
energies depends on the relevant quantum numbers. Moreover, the single-electron wave functions can 
be influenced by the presence of this additional spatially dependent field. We numerically solved the 
Bogoliubov-de Gennes equations with the Hartree-Fock term for a clean metallic nanocylinder and 
found a shift of the curve representing the thickness-dependent oscillations of the critical temperature 
(the energy gap, the order parameter etc.) to larger diameters. Though the difference between the 
superconducting solutions with and without the Hartree-Fock interaction can, for some diameters, 
be very significant, the above mentioned shift is less than typical metallic unit-cell dimensions and, 
so, has no practical worth. This allows one to significantly simplify the problem and, similar to 
bulk, ignore the Hartree-Fock potential when solving the Bogoliubov-de Gennes equations in the 
nano-regime. 

PACS numbers: 74.78.-w, 74.78.Na 



I. INTRODUCTION 

Advances in nanofabrication technology resulted re- 
cently in high-quality metallic superconducting ultrathin 
nanofilmsi^ and nanowires4^^ In most samples the 
electron mean free path was estimated to be about or 
larger than the nanofilm/nanowire thickness ^diL In this 
case the effects of the transverse quantization are not 
shadowed by impurity scattering and, hence, the con- 
duction band splits up into a series of single-electron 
subbands resulting from the quantized transverse modes. 
This will have a pronounced effect on the superconduct- 
ing properties (see, for instance, Refs. H and d and ref- 
erences therein). Notice that high-quality nanofilms do 
not exhibit significant indications of defect- or phase- 
driven suppression of superconductivity (see discussion in 
Ref. [2|). For high-quality nanowires the phase-fluctuation 
effects were shown to seriously influence the supercon- 
ducting state only in narrowest aluminum specimens with 
width w 5 — 8 nm4ii£ Thus, the transverse quantum 
confinement is the major mechanism governing the su- 
perconducting properties in this case. Therefore, it is 
timely to study in a more detail a clean nanoscale super- 
conductor in the presence of quantum confinement. 

Quantum confinement breaks the translational sym- 
metry and, so, the superconducting order parameter be- 
comes position dependent. The well-known BCS ansatz 
for the ground state wave function is not applicable in 
this case, and the Bogoliubov-de Gennes (BdG) equa- 
tions are a relevant tool to investigate equilibrium super- 
conducting properties. Recent numerical studies of the 
BdG equations for nanofilms^ and nanowire s 9 ' 11 ! 12 show 



that the transverse quantum confinement has a substan- 
tial impact on the superconducting solution. However, 
the BdG equations investigated in Refs. d, [ll| and 
12, were solved without the Hartree-Fock (HF) potential. 
The reason is that in bulk, the superconducting solution 
is not sensitive to the HF term in the BdG equations 13 , 
and one can assume that a similar conclusion holds for 
the broken translational symmetry. However, at present 
there is no detailed investigations on this subject and, so, 
such a study is needed. 

In the bulk BdG equations, the HF potential is not 
spatially dependent and, so, it produces the same contri- 
bution to all single-electron energies, with no dependence 
on the relevant quantum numbers. Hence, the Fermi level 
(the chemical potential) acquires the same contribution, 
as well, and the single-electron energies measured from 
the Fermi level are not changed. It is well-known that 
the BdG equations are derived within the grand canoni- 
cal formalism and, so, the electron energies appearing in 
the basic expressions absorb the chemical potential. As a 
result, the superconducting solution is insensitive to the 
HF potential. The situation is different in the presence 
of quantum confinement. The translational symmetry is 
now broken, the HF mean field is position dependent, 
and, so, its contribution to the single-electron energies 
is a function of the relevant quantum numbers. Fur- 
thermore, the single-electron wave functions themselves 
are influenced by the presence of the HF field, i.e., an 
additional spatially-dependent potential. Therefore, one 
can expect that the HF term in the BdG equations can 
change the superconducting solution in the presence of 
quantum confinement. It is of importance to clarify to 
what extent this will be through. In particular, this con- 
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cerns the thickness-dependent oscillations (i.e., quantum- 
size oscillations) of the superconducting properties typi- 
cal of high-quality nanofilms and nanowiresi^^. 

In the present work, based on the particular case of 
a superconducting clean metallic nanowire of the cylin- 
drical form, we compare the superconducting numerical 
solutions of the BdG equations with and without the HF 
potential. We find that these solutions are indeed dif- 
ferent and, for some nanocylinder diameters, this can be 
very significant. However, a more close look reveals that 
this difference is actually expressed in a small shift of the 
curve representing thickness-dependent oscillations of the 
critical temperature (or other important superconduct- 
ing quantities) up to larger diameters. The shift is less 
than typical metallic unit-cell dimensions and, so, can be 
ignored. 

The paper is organized as follows. In Sec. [TT] the for- 
malism of the BdG equations is outlined, with the fo- 
cuser on the features related to the cylindrical confining 
geometry. In addition, the Anderson approximate solu- 
tion to the BdG equations is discussed in this section. 
To check the effect of the HF term on the single-electron 
wave functions, the Anderson solution is constructed by 
assuming that these wave functions do not change in the 
presence of the HF interaction. In Sec. IIIII numerical re- 
sults of the BdG equations with and without the HF term 
are discussed. Based on the Anderson solution, here we 
also investigate the effect of the HF term on the single- 
electron wave functions. 



II. BOGOLIUBOV-DE GENNES EQUATIONS 
AND ANDERSON'S RECIPE 



We focus on the basic superconducting properties of a 
metallic clean cylindrical nanowire (with diameter D — 
2R and length L) in the quantum-size regime when the 
transverse quantization of the single-electron spectrum is 
of importance. In the presence of quantum confinement 
the translational invariance is broken, and the order pa- 
rameter appears to be position-dependent, i.e., A(r). It 
is well-known that the BdG equations are a common and 
useful approach to investigate such a situation. Gener- 
ally, these equations can be represented as follows: 



H & \u v )+&\v v ), 

k*\u u )-H;\v v ), 
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(lb) 



where E v stands for the quasiparticle energy, \u u ) and 
\v v ) are the particle-like and hole-like ket vectors. In the 
clean limit the single-electron Hamiltonian in Eqs. (TTa|) 
and (llb[) is of the form [for zero magnetic field, A = 0] 



h p = m 



2m, 



+ $HF(i) + V conl (r)- E F , (2) 



with f and p the position and momentum operators, Ep 
the Fermi level, m e the electron band mass (set to the free 
electron mass), V^ on f(r) the confining interaction, and 



&hf(y) the HF potential. In bulk the confining inter- 
action can be neglected and we arrive at the usual BCS 
picture based on plane waves. Below we adopt the sim- 
plest choice of the confining interaction potential: zero 
inside and infinite outside the wire. The gap-operator A 
in Eqs. ([Taj) and (JTbJ) is related to the order parameter 
by A = A(f). 

As a mean-field theory, the BdG equations are solved 
in a self-consistent manner with the self-consistency re- 
lations given by 

A(r)=gJ2^\uu)(vu\T)[l-2f v ], (3a) 
i>ec 

- -g]T[KiK)| 2 /, + |(rK)| 2 (l - /„)] , (3b) 

where g > is the coupling constant, f v = l/(ef E " + 1) 
is the Fermi function [j3 — 1/ (fcgT) with T the tempera- 
ture and fee the Boltzmann constant]. In Eq. (|3a[) C in- 
dicates the set of quantum numbers corresponding to the 
single-electron energy ^ v (measured from the Fermi level) 
located in the Debye window £„ e c <= h^c, fu^D\ (wd is 
the Debye frequency) , where f v absorbs the HF potential, 
i.e., 



i„ = {u v \H e \u v 



mh:\ Vv ) 



(4) 



The cut-off in Eq. (|3"a|) is known^ to be a payment for 
using a simplified delta-function approximation for the 
electron-electron interaction. Such a regularization is not 
needed in Eq. (|3b[) . For our confining interaction (i.e., 
zero inside and infinite outside) we have 

(r\u v ) ={t\v v ) =0 (5) 

at the sample surface, i.e., r <E S. Periodic boundary 
conditions with unit cell L can be applied in the direction 
parallel to the nanowire. 

The Fermi level (i.e., the chemical potential) is deter- 
mined from 

2 r 

n e = 2 ^2\{u v \u v )f v + (v u \v v )(l - /„ 



(6) 



where n e is the mean electron density. We use the BdG 
equations in the parabolic band approximation and, so, 
as discussed in Ref. [H, an effective Fermi level should be 
introduced, to recover the correct period of the quantum- 
size oscillations. For aluminum (the aluminum parame- 
ters are used below) Ep = 0.9 eV for D w 10 nm (see 
Ref.@). For D ~ 1 — 2 nm, Ep shifts systematically from 
this value up, due to Eq. ©. 

Due to the chosen confining geometry, it is convenient 
to use cylindrical coordinates p, tp and z. In this case the 
order parameter (the anomalous pairing potential) and 
HF mean field (the normal potential) depend only on the 
transverse coordinate, i.e., A(p) and &hf(p), and (r\u„} 
and (r|i;„) are represented in the form (y = {j,m,k}) 



{r\Uj m k) 
(r\Vjmk) 



2^ VI 



Ujmk(p) 
Vjmk(p) 



(7) 
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with j controlling the number of nodes in the transverse 
direction, m the azimuthal quantum number, and k the 
wave vector of the quasi-free electron motion along the 
nanocylinder. Inserting Eq. (0 into Eqs. (fla|) and (|lb[) . 
we recast the BdG equations as 

[Ej m k ~£p- ®hf(p)] Uj m k(p) = A{p)v 3m k(p), (8a) 
[Ejmk + C p + <5>Hf(p)] Vjmk(p) = A(p)ujmk(p), (8b) 

where A(p) is real, and 



£o = — 



-dp 2 



ld_ 

pdp 



- k 2 )-E F . 



(9) 



2m e v op* p op p 
The self-consistency relations can be rewritten as, 

u jmk {p)v 3 mk (p) [l - 2fjmk] , (10a) 



jmkEO 



<S>hf(p)- 
9 v^l 



jmk 



[ujmMfjmk + V 2 jmk {p)(l - fj mk ) , (10b) 



Notice that it is possible to deal with Anderson's recipe, 
invoking the eigenfunctions of C p + $hf(p)- However, 
below we are interested in Eq. (TH| because it helps to 
clarify how a change in the single-electron wave func- 
tions due to the HF potential, can contribute to the 
problem of interest. To be accurate, the Anderson ap- 
proximation should be based on the true single-electron 
wave functions. We recently found that in this case 
the error in Anderson's solution for D < 2 — 3 run is 
less than one-two percentsJ^ Hence, comparing the re- 
sults of numerically solving Eqs. (|8a|) and (|8b|) with 
the data based on Eq. (fTi|) . we can reach unambigu- 
ous conclusions about the role of the changes in the 
single-electron wave functions due to the HF interac- 
tion. As follows from Eq. (TTJJ [see, for instance, Ref . [l6| . 
the Anderson-approximation results in the BCS-like self- 
consistent equation 

9fm', jm A ^ [l _ 2f]mk ^ (15) 



j mh 



kec J L mk + A,- 



with 



with Uj m k{p) and Vj m k(p) real. To numerically solve 
Eqs. (jSaj) and (|8bj) . we expand the transverse particle- 
like and hole-like wave functions as 



Ujmk(p) \ _ ST^ 
V jm k{p) J 



Vjm(p), (11) 



with 



^Jm(p) = ^7 rJ m (aj m ^), 

iiJm+l\ a Jm) H 



(12) 



where J m {x) is the Bessel function of the first kind of 
the m-order, and aj m is the Jth zero of this function. 
This allows one to convert Eqs. ([8a| and (|8bj) into a ma- 
trix form. Then, a numerical solution can be obtained 
by diagonalizing the corresponding matrix, and self- 
consistency is reached by iterating Eqs. (jlOap and (|10b[) . 
One should keep in mind tha"bii (u„|u„) + (v u \v u ) = 1 
and, so, 



d PPi u %nk(p) + Vj mk {p)] 



(13) 



As seen, the vector (uj m k,j;Vjmk,j)' T {J = 0,1,...) is 
normalized. 

In addition to the above procedure, below we use 
the Anderson approximate solution, as welli^. Within 
this approximation, instead of the expansion given by 
Eq. (fTTj) . it is assumed that 

Uj m k(p) = Uj mk 'dj m (p), Vj m k{p) = Vjmkftjm{p)- (14) 



Equation (I14|) means that we seek a minimum of the BdG 
thermodynamic functional in the subspace of Uj m k{p) 
and Vj m k(p) proportional to the eigenfunctions of C p . 



and the interaction-matrix element given by 
R 

_ a 

Qj'm',jm — 



2ttL 







£,jmk 



2m e 



R 2 



where 
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Inserting Eq. (|3b|) into Eq. (fl9|) . one obtains 



jmk 



Qj'm' jm 



l - 



Cjmk (1 ^fjmk^ 



jmk 



A 2 



(16) 



(17) 



For the single-electron energy appearing in Eq. (|15[) we 
have 

v 2 



(18) 



(19) 



(20) 



We should not forget about Ep appearing in the single- 
electron energy given by Eq. (fT8|) . It is fixed through 
Eq. ([6]) that is now of the form 



ttR 2 L E I 1 



^jmk (1 jmk ) 



jmk 



jmk 



A 2 



(21) 



Thus, in the Anderson approximation introduced by 
Eq. (fTi|), one needs to solve Eqs. fT5]) and (|2"0|) . keep- 
ing Eq. (|2ip . As already mentioned above, comparing a 
numerical solution of Eqs. (|8"aj) and (|8b[) with the solution 
based on Anderson's recipe, we can check the effect of the 
HF interaction on the single-electron wave functions. 
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D (jam) D (nm) 

FIG. 1: (a) The spatially averaged superconducting order parameter A/Ab u ik and (b) critical temperature T c /T Ci buik versus 
the nanowires diameter D as calculated from the BdG equations (I8a|l and (|8bj) at zero temperature. 



III. NUMERICAL RESULTS 

In this section we investigate and discuss numerical 
self-consistent solutions of Eqs. (|5a|) and (|8b]) with the 
HF potential (the full version) and without it (the trun- 
cated version, by setting &hf(p) = in the relevant 
expressions). Results are also compared with a solution 
of Eqs. (Unj) and (f2T))) . All the calculations are performed 
with the parameters typical for aluminu m 13 ' 14 : fiioo = 
32.31 meV; gN(0) = 0.18, with N(0) = m e k F /(2n 2 h 2 ) 
the bulk density of single-electron states at the Fermi 
level [for Ep see discussion after Eq. ©] and k F the bulk 
Fermi wavevector. For these parameters the bulk BCS 
coherence length £ = 1.6 is significantly larger than 
the nanocylinder diameter. However, contrary to the 
ordinary Ginzburg-Landau picture, the superconducting 
order parameter now exhibits significant spatial varia- 
tions in the transverse direction due to the broken trans- 
lational symmetry. The length of the nanocylinder is 
taken as L = 1 ^m> Xp = 2ir/kp. This is an optimal 
choice, upholding, on one side, the use of periodic bound- 
ary conditions in the z direction and, on the other side, it 
results in a reasonable calculational time. As opposed to 
the truncated BdG equations, their full version requires 
much more time for convergence of the numerical proce- 
dure, and this time increases proportionally with L 2 . Nu- 
merically solving the Anderson equations (fT5|) and (j2"0|) 
is less time-consuming and, so, we take L — 5 /im in this 
case. 

In Fig. [IJa) the spatially averaged order parameter 

R 




o 



calculated from Eqs. (|8a|) and (|8b|) . is plotted in units of 
the bulk order parameter (Ah u ik = 0.25 meV) versus the 
nanocylinder diameter with and without the HF mean 
field. In Fig. [TJb) the corresponding critical tempera- 



ture T c (in units of the bulk one) is given. As seen, 
both data-sets exhibit pronounced size-dependent oscilla- 
tions, typical of high-quality superconducting nanofilms 
and nanowires with uniform thickness 1 ^^^. Such os- 
cillations result from single-electron subbands forming 
due to the transverse quantization of the electron mo- 
tion. With an increase in the nanowire diameter, the 
subbands shift down in energy. Each time when a new 
subband comes into the Debye window around the Fermi 
level, the number of single-electron states contributing 
to the superconducting order parameter increases, and a 
size-dependent superconducting resonance develops. As 
follows from Fig. [TJ the quantum-size oscillations corre- 
sponding to the full version of the BdG equations are 
somewhat shifted up. Hence, the difference between the 
two sets of data is most significant for those diameters, 
where a size-dependent superconducting resonance in the 
case without the HF interaction is already present while 
in the full version such a resonance only starts to develop. 
The difference is not so significant but still survive when 
the resonance comes into its decay stage. When the res- 
onance is fully decayed (the off-resonant regime), the HF 
corrections are practically negligible, and we arrive at the 
situation similar to bulk. Notice that small differences 
between the numerical results of the full and truncated 
BdG equations in the off-resonant regime (due to beat- 
ing patterns of the corresponding curves), are because 
of the chosen nanowire length. Indeed, as follows from 
calculations for several selected off-resonant diameters, 
such beating patterns disappear when L increases up to 
20 — 30 /im, and the results with and without the HF 
interaction approach each other. 

Notice that maxima of T c /T C: buik in Fig. [ljb) are gen- 
erally higher than those of A/A bu ik in Fig. [TJa). This 
is due to formation of new Andreev-type states induced 
by the transverse quantum confinement (see details in 
Ref. [13), which results in a decrease of A/(fcsT c ) below 
the bulk value 1.763 at the resonant points. As seen from 
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FIG. 2: Upper panel: (a) the superconducting order parameter A(p), (b) the local density n e (p) and (c) the distribution 
u ^mk(p) + t 'jm7c(p) versus p/R for D = 2.24 nm. The lower panel: the same but for D = 2.6 nm. The results for the truncated 
(no HF) and full (HF) BdG equations are plotted. 



Fig. [T] inclusion of the HF interaction can slightly reduce 
the resonant enhancements, with practically no effect on 
the ratio A/(k B T c ). 

In Fig.[5]we present different quantities calculated with 
the full and truncated versions of Eqs. (jSaj) and (|8b|) for 
two diameters: the upper panel, for D — 2.24 nm; and the 
lower panel, for D = 2.6 nm. The upper panel represents 
the situation when a superconducting resonance is devel- 
oped for the truncated version but is not yet present for 
the full version of the BdG equations. In Fig. [2ja) the 
superconducting order parameter A(p) calculated with 
(HF) and without the HF interaction (no HF, the in- 
set), is plotted versus the transverse coordinate p/R for 
T = 0. The spatial distribution of the pair condensate is 
very different for these two cases: the data without the 
HF interaction are larger by an order of magnitude, and 
even the profile of A(p) is different. In Fig. ©(b) the 
local electron density, i.e., 

n e{p) = ^ [u 2 Jmk {p)fjmk + Vj mk (p)(l - f Jm k) , 

jmk 

is shown for the same diameter. As can be expected, 



now the difference between the two data-sets is not so 
significant (we keep to the same value of the mean elec- 
tron density n e ). Due to the attractive character of the 
effective electron-electron interaction, the HF potential 
forces electrons to go closer to the nanocylinder cen- 
ter. However, the confining interaction has the major 
effect on n e (p) as compared to the HF potential pro- 
ducing only some small corrections. From the results 
for the local electron density, it is possible to expect 
that the single-electron wave functions are also not very 
sensitive to the HF interaction. For nanowires, |iti/(r)| 
and |iv(r)| is nearly proportional to the corresponding 
single-electron wave function [see discussion above, af- 
ter Eq. ((HI)]. Hence, due to Eq. (fP3]l . the quantity 
u "jmk(p) + v jmk(p) can provide us with the informa- 
tion about the single-electron distribution. In Fig. [2fc) 
u "jmk(p) + v jmk(p) ^ s plotted versus p/R for the quantum 
numbers most sensitive to including the HF interaction. 
We can indeed see that the effect of the HF potential on 
the wave functions is minor. Similar conclusions can be 
obtained from the lower panel of Fig. O The only ex- 
ception is that the superconducting order parameter in 
Fig. [2]Jd) [D = 2.6 nm] does not change so much when 
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FIG. 3: A/Abuik versus the nanowire diameter (T = 0): tri- 
angles correspond to the Anderson approximation [the HF 
field is included]; stars and squares are the numerical results 
of the full and truncated BdG equations, respectively. 
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FIG. 4: HF versus no HF: A/A bu ik as a function of the 
nanowire diameter D at zero temperature [squares are the 
results of the truncated BdG equations (no HF interaction); 
triangles are the results of the Anderson approximation, see 
Eqs. (HU and dl). 

including the HF potential. Notice that n e (p) given in 
Fig. [^e) is practically the same as in Fig. HJb). How- 
ever, this is not true for u 2 mk (p) + Vj mk (p) [compare 
panel (c) with panel (f)]. The point is that the inte- 
gral J Q dppn e (p) — n e R 2 /2 changes with the radius but 
for the single-electron distribution u 2 mk (p) + v 2 mk (p) we 
have Eq. (p])]. 

From the results presented in Fig. [21 one expects minor 
effects on the single-electron wave functions due to the 
incorporation of the HF interaction . This expectation 
can be put on a more solid ground by using the An- 
derson approximation based on Eqs. (fTS")) and ([!"§)) . We 
remind that the Anderson approximation is quite good 
for superconducting nanowires provided that it involves 



the true single-electron wave functions. Equations (|15p 
and ([IT))) follow from Eq. (fT4"|) and, hence, as assumed, 
the single-electron wave functions are not altered by our 
position-dependent HF interaction. If this is a reason- 
able assumption, results of the Anderson approximation 
constructed in this way, should be close to the results 
of the full BdG equations. As seen from Fig. [31 this is 
indeed the case. We can conclude that the thickness- 
dependent shift of the superconducting resonances in the 
presence of the HF interaction has nothing to do with the 
single-electron wave functions. Its mechanism is due to 
the fact that the position-dependent HF potential results 
in a change of the single-electron energies measured from 
the Fermi level. 

So far we considered extremely narrow nanowires, for 
the sake of simplicity. However, a similar shift (w 0.01 — 
0.02 nm) of the quantum-size oscillations due to the HF 
term survives until the total decay of the quantum-size 
oscillations (up to diameters of about 50 — 70 nm) . In 
particular, such a shift is clearly seen in Fig. [4l where 
numerical results of the truncated BdG equations for 
D = 10 — 10.5 nm are compared with a solution of 
Eqs. (fl"5j) and([Tn]) including the HF potential. Thus, we 
arrive at the following picture. In the vicinity of a super- 
conducting resonance, the bottom of some single-electron 
subband is situated close to the Fermi level. Therefore, a 
repositioning of this subband with respect to the Fermi 
level can result in a significant change of the number of 
single-electron states in the Debye window and, so, in a 
remarkable increase/decrease of superconducting charac- 
teristics. However, when bottoms of all single-electron 
subbands are quite apart from the Fermi level, i.e., in 
the off-resonant regime, a move of these subbands in en- 
ergy produces much less important effect on the number 
of single-electron states located in the Debye window. 
This is why the decay of a superconducting resonance is 
accompanied by a depletion of the influence of the HF 
potential. 



IV. CONCLUSIONS 

Quantum confinement breaks the translational sym- 
metry in nanostructured superconductors. In this case, 
despite the delta-function approximation for the electron- 
electron interaction, the HF potential becomes position 
dependent, and its contribution to the single-electron en- 
ergy (measured from the Fermi level) is a function of the 
relevant quantum numbers (contrary to bulk!). By nu- 
merically solving the Bogoliubov-de Gennes equations for 
a clean metallic nanocylinder, we have shown that such 
a feature results in a shift of the curve representing the 
thickness-dependent oscillations of the critical tempera- 
ture (the energy gap, the order parameter etc.) to larger 
diameters. Though difference between the numerical re- 
sults with and without the HF potential can be rather 
significant for some given diameters, the above mentioned 
shift is less than typical metallic unit-cell dimensions and, 
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so, can be ignored. In particular, smoothing due to the 
thickness fluctuations in real samples will shadow this 
effect. Thus, when solving the Bogoliubov-de Gennes 
equations in the nano-regime, one is able to neglect the 
HF term, which significantly simplifies a numerical pro- 
cedure. 
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